數值微分標誌著從微積分無限平滑的理論世界，轉向離散且有限的數位運算環境。我們以可量測的步長 $h$ 取代無限小的極限。雖然函數 $f$ 在 $x_0$ 處的理論導數定義為 $$f'(x_0) = \lim_{h \to 0} \frac{f(x_0 + h) - f(x_0)}{h}$$，但電腦系統無法直接計算極限。相反地，我們使用有限差分公式，並產生一個可量化的代價，稱為 截斷誤差。
1. 导数的幾何意義
為了近似 $f'(x_0)$，我們觀察鄰近的點。根據方向選擇的不同，可推導出兩種主要公式：
- 前向差分公式： 當 $h > 0$ 時使用。它向前觀察 $x_0 + h$。
- 後向差分公式： 當 $h < 0$ 時使用。它向後觀察 $x_0 + h$（此時 $h$ 為負數）。
在現實工程應用中，例如計算曲線軌跡的弧長時，我們經常依賴這些近似方法：$$L = \int_{0}^{48} \sqrt{1 + (f'(x))^2} dx = \int_{0}^{48} \sqrt{1 + (\cos x)^2} dx$$ 若 $f(x)$ 僅在離散的感測器點上已知，則數值微分是唯一可行的途徑。
2. 透過插值的數學推導
為了近似 $f'(x_0)$，首先假設 $x_0 \in (a, b)$，其中 $f \in C^2[a, b]$，且 $x_1 = x_0 + h$。我們構造由 $x_0$ 與 $x_1$ 決定的第一個拉格朗日多項式 $P_{0,1}(x)$：
步驟一：插值函數建構
$f(x) = P_{0,1}(x) + \frac{(x - x_0)(x - x_1)}{2!} f''(\xi(x))$
步驟二：求導
對兩側求導並於 $x = x_0$ 評估，可得基本關係式：
$$f'(x_0) = \frac{f(x_0 + h) - f(x_0)}{h} - \frac{h}{2} f''(\xi)$$
3. 誤差項與收斂性
項 $-\frac{h}{2} f''(\xi)$ 即為我們的截斷誤差。此公式顯示準確度為 $O(h)$，表示若將步長 $h$ 減半，誤差也大致減半。然而我們必須謹慎：雖然較小的 $h$ 可減少截斷誤差，但最終會導致 捨入誤差 因分子中近乎相等的數字相減所導致。
🎯 核心原則：有限差分
數值微分以有限弦線取代極限。我們近似的品質完全取決於步長 $h$ 與函數的平滑程度（即二階導數）。
$f'(x_0) \approx \frac{f(x_0+h) - f(x_0)}{h}$，誤差界為 $\frac{h}{2} \max|f''(\xi)|$